Condensation of N interacting bosons: Hybrid approach to condensate fluctuations 
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We present a new method of calculating the distribution function and fluctuations for a Bose- 
Einstein condensate (BEC) of N interacting atoms. The present formulation combines our previous 
master equation and canonical ensemble quasiparticle techniques. It is applicable both for ideal and 
interacting Bogoliubov BEC and yields remarkable accuracy at all temperatures. For the interacting 
gas of 200 bosons in a box we plot the temperature dependence of the first four central moments 
of the condensate particle number and compare the results with the ideal gas. For the interacting 
mesoscopic BEC, as with the ideal gas, we find a smooth transition for the condensate particle 
number as we pass through the critical temperature. 



Bose-Einstein condensation (BEC) near the critical 
temperature T c is a fascinating subject. Indeed the time 
line from the early studies of Uhlenbeck p| concerning 
the unphysical cusp at T c in Einstein's analysis, to the 
"Grand canonical catastrophe" 0], arising from ill fated 
attempts to describe fluctuations [3j, covers most of the 
twentieth century. Naturally, a major stimulus for the 
present work is the pioneering BEC experiments in liq- 
uid He-II and in laser cooled gases pj- The attendant 
theoretical advances in the study of BEC in ideal and 
interacting gases likewise constitute a rich field 0, • 

Furthermore, subtle issues such as the statistics of con- 
densate atoms in a mesoscopic system of N ~ 10 2 — 10 3 
particles are now of interest. Condensate fluctuations can 
be measured by means of a scattering of series of short 
laser pulses 0, see also 0. We note that the BEC is 
often referred to as an atom laser; indeed the problem 
of BEC statistics near T c is analogous to studying the 
photon statistics of the laser in the passage from below 
to above threshold 0, 0] . 

In the present Letter, we give, for the first time, a 
simple, surprisingly accurate account of fluctuations in 
an interacting Bogoliubov gas valid for all temperatures. 
We omit, however, effects of interaction between Bogoli- 
ubov quasiparticles since we treat a weakly interacting 
gas. The analysis is based on a master equation approach 
deriving from a union of the ideal gas (in the spirit of the 
quantum theory of the laser) BEC treatment of the first 
two papers in this series on the condensation of N bosons 
[To| and the results for the interacting gas BEC obtained 
based on the canonical ensemble quasiparticle formalism 



We emphasize that although useful papers have been 
published dealing with various limiting cases, so far there 
has been no treatment of this problem valid at all tem- 
peratures 0|- Ref. 0] presented analytical formulas for 
all moments of the condensate particle number fluctua- 
tions in the weakly interacting Bose gas. However, 0] 
and other approaches (e.g. Ref. 0]) are only valid pro- 
vided the average number of condensate particles is much 
larger than its variance. However, near and above T c this 
is not true, and this causes the failure of such treatments. 

For an ideal Bose gas the master equation approach 



of [lfj naturally includes the N particle constraint and 
provides an analytical solution for the partition function 
and fluctuations accurate at all temperatures, as shown 
in Fig. ^ F° r an interacting gas, the problem of fluctu- 
ations is rather delicate. In a recent paper 1 1 ol | we gave 
a preliminary master equation analysis for N interacting 
atoms which accurately described the average number 
of particles in the condensate. However the fluctuations 
were handled less well. The present analysis gives all 
central moments with remarkable accuracy when com- 
pared to 0| below T c . Above T c the present results go 
smoothly into the ideal gas limit as they must. 




FIG. 1: The variance for the condensate particle number as a 
function of temperature for an ideal gas of N = 200 particles 
in an isotropic harmonic trap obtained by the master equation 
approach of 10] (solid line), "exact" numerical simulations in 
the canonical ensemble (dots) and the grand canonical answer 
(dashed line). Result of Refs. Il2| is plotted as dash-dot 
line. Small dots show the thermodynamic limit formula of 
Politzer Q. 

The central tool used in the ideal gas analysis of 0| 
and the interacting gas study of was the laser-like 
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master equation for the probability p no of finding no 
atoms in the condensate, given that there are N total 
particles 

1 

-Pn = -K nQ (n + l)Pn + K nQ _ 1 n p no - 1 - 



p n = n"e~"/n!, as would be expected for a coherent 
state. 

Using the distribution function © we find that, in the 
validity range of [12 (at low T), the first thr ee central 
moments fj, m =< (no — no) m > are 

n = N-H, u 2 = (1 + T))H + aH 2 , (7) 



H no n p no + H no+1 {n + l)p„ 0+ i, (1) 

where n is an uninteresting rate constant, H no and K no 
are heating and cooling coefficients. In equilibrium the 
rates of any two opposite processes are equal to each 
other, e.g., K no (n + l)p m = H na+1 (n + l)p no+1 . The 
detailed balance condition yields 

Pn + 1 _ Kng ^ 
Pno H na 4. 1 

Since the occupation number of the ground state can- 
not be larger than TV there is a canonical ensemble con- 
straint pn+i = and, hence, Kn — 0. In contrast to 
p no , the ratio p no +i/p no as a function of n shows simple 
monotonic behavior. We approximate K no and H no by a 
few terms of the Taylor expansion near the point no = N 

K no = (N -n )(l + n)+a(N -n ) 2 , (3) 



H no =H+(N-n ) V + a(N-n ) 2 . (4) 

Parameters H., r\ and a are independent of no; they are 
functions of the occupation of the excited levels. We 
derive them below by matching the first three central 
moments in the low temperature limit with the result of 
[T^ . We note that the detailed balance equation J2J is the 
Pade approximation of the function p no+ i/p no . Pade 
summation has proven to be useful in many applications, 
including condensed-matter problems and quantum field 
theory. 

Eqs. yield an analytical expression for the con- 

densate distribution function 

= I (iV-no-l + a;i)!(JV-no-l + a; 2 )! 
Pn ° Z N (N-n )\(N-n + (l + r])/a)\ ' U 



where xi t 2 = (tj ± \Zv 2 4a7i)/2a and Zn is the nor- 

N 

malization constant determined by ^ p„ — 1. In the 

n =0 

particular case 77 = a = Eq. J5j reduces to 

■y 'j.fN—na 
Pno — 1y 7T7 sTj (6) 

Z N (N - n )! 

where Zn = e n T(N + 1,H)/Nl is the partition function 
and r is an incomplete gamma-function. For an ideal gas 
Eq. © describes accurately the condensate statistics at 
low temperature [ljj. The statistics is not Poissonian 



"3 = -H(l + r ] + aH)(l + 2r 1 + 4aH). (8) 
Eqs. 0, © thus yield 

„_J(s_, + te), (9) 

-sG-w-e)- (10) 

On the other hand, the result of 01 f° r an interacting 
Bogoliubov gas is (see also 0] for no and n 2 ) 

n = N-J2[(ul + v 2 k )f k + vl], (11) 
M2 = £ [(1 + 8u 2 v 2 )(f 2 + f k ) + 2u 2 v 2 ] , (12) 

fc^O 

Ma = - E K + [C 1 + 16u fe^)( 2 /fe + 3 A 2 + /*)+ 

fe#0 

4^(1 + 2/ fe )], (13) 

where = l/[exp(Ek/kBT) — 1] is the number of ele- 
mentary excitations with energy E k present in the sys- 
tem at thermal equilibrium, Uk and v k are Bogoliubov 
amplitudes. Substitute for fio, un and 113 in Eqs. JSJ, 
(fTUfl their expressions of Ref. [12] lll|) - (|13|l yields the 
unknown parameters TL, r\ and a. The beauty of our 
"matched asymptote" derivation is that the formulas for 
TL, 77 and a are applicable at all temperatures, i.e. not 
only in the validity range of [12j. The distribution func- 
tion J5| together with Eqs. JDJ, I jlOjl provides complete 
knowledge of the condensate statistics at all T. Taking 
Vk = and u k — 1 in l|ll|) - (|13|l we obtain the ideal gas 
limit. 

Next we test our method for an ideal gas (in a har- 
monic trap). In this case an "exact" numerical simula- 
tion in the canonical ensemble is available for com- 
parison. Results of such simulation are shown by dots in 
Figs. [21 and In Fig. |S]we plot the distribution func- 
tion for the number of atoms in condensate at different 
temperatures and N = 200. At T < T c the distribu- 
tion shows a sharp peak near no and becomes broader at 
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higher T. The present Eq. (JSJ) (solid line) yields excel- 
lent agreement with the "exact" dots at all temperatures. 
Fig. shows the average condensate particle number no, 
its variance, third and fourth central moments /i m and 
fourth cumulant K4 Jl7| as a function of T for N = 200 
particles in a harmonic trap. Solid lines are the result of 
the present approach (we call it CNB5 0|) which is in 
remarkable agreement with the "exact" dots at all tem- 
peratures both for ji m and K4. Central moments and 
cumulants higher than fourth order are not shown here, 
but they are also remarkably accurate at all tempera- 
tures. Results of are given by dashed lines which are 
accurate only at sufficiently low T . Deviation of higher 
order cumulants (to = 3,4,...) from zero indicates that 
the fluctuations are not Gaussian. 



N=200 




FIG. 2: Distribution function for the condensate particle 
number no at different temperature obtained by the present 
approach (solid lines) and "exact" numerical simulations in 
the canonical ensemble (dots). The results are obtained for 
an ideal gas of N = 200 particles in an isotropic harmonic 
trap. 

Clearly our method passes the ideal gas test with fly- 
ing colors. Please note the excellent agreement with the 



exact analysis for the third central moment and fourth 
cumulant K4 given in Fig. |31 

Next we apply the present technique to N interact- 
ing Bogoliubov particles confined in a box of volume V. 
The interactions are characterized by the gas parame- 
ter an 1 / 3 , where a is the s-wave scattering length and 
n = N/V is the particle density. The energy of Bo- 
goliubov quasiparticles Ek depends on no, hence, the 

JV 

equation no — n oPn for no must be solved self- 

no=0 

consistently. In Fig. 0] we plot no, the variance Ano, 
third and fourth central moments as a function of T for 
an ideal and interacting (an 1 / 3 = 0.1) gas in the box. 
Solid lines show the result of the present approach, while 
[l2T | is represented by dashed lines. The present results 
agree well for all u™. with |l2j| in the range of its validity. 
Near and above T c [12] becomes inaccurate. However, the 
results of the present method are expected to be accurate 
at all T. Indeed, in the limit T ^> T c the present results 
(unlike [T^]) merge with those for the ideal gas. This is 
physically appealing since at high T the kinetic energy 
becomes much larger than the interaction energy and the 
gas behaves ideally. Similar to the ideal gas, the interact- 
ing mesoscopic BEC no(T) exhibits a smooth transition 
when passing through T c . 

One can see from Fig. 0] that the repulsive interac- 
tion stimulates BEC, and yields an increase in no at in- 
termediate temperatures, as compared to the ideal gas. 
This effect is known as "attraction in momentum space" 
and occurs for energetic reasons [l^. Bosons in differ- 
ent states interact more strongly than when they are in 
the same state, and this favors multiple occupation of a 
single one-particle state. 

In conclusion, in this paper we presented a simple 
method which, for the first time, yields an accurate de- 
scription of the distribution function and fluctuations for 
mesoscopic interacting Bogoliubov BEC in the canonical 
ensemble at all temperatures. Our approach combines 
the analytical results of [12j with the laser-like master 
equation of 
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FIG. 3: Average condensate particle number < no >, its variance Ano = y < ( n o — no) 2 >, third and fourth central moments 
< (no — no) m > (m = 3,4) and fourth cumulant K4 as a function of temperature for an ideal gas of N = 200 particles in a 
harmonic trap. Solid lines (CNB5) show the result of the present approach. [l2Tl yields dashed lines (CNB3). Dots are "exact" 
numerical simulation in the canonical ensemble. The temperature is normalized by the thermodynamic critical temperature 
for the trap T c = HljN 1/3 /k B ((3) 1/3 , where u> is the trap frequency. 
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FIG. 4: Average condensate particle number, its variance, third and fourth central moments as a function of temperature for 
an ideal (on 1//3 = 0) and interacting (an 1 / 3 = 0.1) Bose gas of N = 200 particles in a box. Solid lines are the result of the 
present approach. |l2| yields dashed lines (CNB3). The temperature is normalized by the thermodynamic critical temperature 
for the box T c = 27rft 2 n 2/3 /fc s MC(3/2) 2/3 , where M is the particle mass. 



